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ABSTRACT 

This is the third in a series of papers studying the astrophysics and cosmology of 
massive, dynamically relaxed galaxy clusters. Our sample comprises 40 clusters iden¬ 
tified as being dynamically relaxed and hot (i.e., massive) in Papers I and II of this 
series. Here we consider the thermodynamics of the intracluster medium, in particu¬ 
lar the profiles of density, temperature and related quantities, as well as integrated 
measurements of gas mass, average temperature, total luminosity and center-excluded 
luminosity. We fit power-law scaling relations of each of these quantities as a func¬ 
tion of redshift and cluster mass, which can be measured precisely and with minimal 
bias for these relaxed clusters using hydrostatic arguments. For the thermodynamic 
profiles, we jointly model the density and temperature and their intrinsic scatter as a 
function of radius, thus also capturing the behavior of the gas pressure and entropy. For 
the integrated quantities, we also jointly fit a multidimensional intrinsic covariance. 
Our results reinforce the view that simple hydrodynamical models provide a good 
description of relaxed clusters outside their centers, but that additional heating and 
cooling processes are important in the inner regions (radii r < 0.5r25oo ~ 0.15 rsoo). 
The thermodynamic profiles remain regular, with small intrinsic scatter, down to the 
smallest radii where deprojection is straightforward (~ 20kpc); within this radius, 
even the most relaxed systems show clear departures from spherical symmetry. Our 
results suggest that heating and cooling are continuously regulated in a tight feedback 
loop, allowing the cluster atmosphere to remain stratified on these scales. 
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1 INTRODUCTION 

Highly dynamically relaxed galaxy clusters represent a mi¬ 
nority of the cluster population, but provide a key laboratory 
for learning about the physics of the intracluster medium 
(ICM) and its interaction with cluster galaxies and active 
galactic nuclei (AGN). Uniquely in these morphologically 
regular systems, the three dimensional, radial profiles of the 
thermodynamic properties of the ICM can be reconstructed 
with minimal uncertainties from projection effects. In addi¬ 
tion, X-ray data can provide precise constraints on the mass 
profiles of these clusters under the assumption of hydrostatic 
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equilibrium. The latter feature comes at the expense of in¬ 
troducing some systematic uncertainty due to non-thermal 
support, but both simulations and direct calibration using 
weak gravitational leasing data show that the overall bias in 
Chandra X-ray mass estimates for relaxed clusters is small 
( <, 10 per cent; Nagai et al. 2007; Applegate et al. 2016). 

Constraints on mass profiles, in comparison with esti¬ 
mates of the mass within a single characteristic radius, en¬ 
able a wider range of investigations. In general, they provide 
a way to identify dynamically comparable radii across clus¬ 
ters, yielding the most meaningful comparison of interesting 
features. More specifically, they allow the scaling relations of 
thermodynamic quantities to be investigated as a function of 
radius. Violations of the self-similar scaling predicted from 
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purely gravitational spherical collapse shed light on the as¬ 
trophysics of the ICM, in particular the heating mechanism 
that prevents classical cooling flows from forming in relaxed 
clusters. While numerous studies have investigated the scal¬ 
ing relations of global cluster properties (for a review, see 
Giodini et al. 2013), there has not yet been a complete study 
of the scaling of cluster thermodynamic profiles, including 
the radial structure of their intrinsic scatter (though see Say¬ 
ers et al. 2013). 

Other papers in this series have focused on identifying 
a sample of massive, relaxed clusters (Mantz et al. 2015), 
exploiting measurements of their gas mass fractions at inter¬ 
mediate radii to constrain cosmology (Mantz et al. 2014b), 
and measuring the average bias of X-ray hydrostatic masses 
for these clusters (Applegate et al. 2016; hereafter, respec¬ 
tively, Papers I, II and IV). Here we assume a concordance 
ACDM cosmology (with parameters h = 0.7, fim = 0.3 and 
Ha = 0.7) and present the thermodynamic profiles of our 
relaxed cluster sample, their scaling properties, and some re¬ 
lated astrophysical results. We briefly review aspects of our 
sample selection and spectral analysis methods in Section 2. 
Section 3 discusses the additional statistical methodology 
used in this work, in particular for performing multivari¬ 
ate regression with intrinsic scatter and/or arbitrary mea¬ 
surement error distributions. In Section 4, we present the 
thermodynamic profiles of the cluster sample and discuss 
some of the immediate consequences of these measurements 
that are available without further model fitting. Section 5 
presents scaling relations of the profiles, including the radial 
structure of their intrinsic scatter, as well as “traditional” 
scaling relations of globally measured quantities. We sum¬ 
marize these results in Section 6. 

In this paper, we follow the convention of defining char¬ 
acteristic cluster masses and radii in terms of the critical 
density of the Universe at the cluster’s redshift, 

Ma = -■nAp„{z)r%, (1) 

where A is often referred to as the “overdensity”. With 

this definition, the self-similar predictions for the scaling of 
the thermodynamic ICM properties can be written (Kaiser 
1986) 

n(rA) oc E{z)^, (2) 

fcT(rA) oc [A^/^£;(^)MA]"^^ 

where n and kT are the number density and temperature 
of the ICM. Here E{z) = H{z)/H(, is the Hubble parameter 

^oc y/ Pct{z)^ , normalized to its value at z = 0 for conve¬ 
nience. Note that, using A as a surrogate for cluster radius, 
the shape of the self-similar density profile is completely 
specified by the expression above, whereas the shape of the 
temperature profile is determined by the shape of the un¬ 
derlying mass profile. Expressions for the self-similar scaling 
of pressure, P = nkT, and pseudo-entropy, K = kTjr?^^, 
follow directly from those above, as do the predicted scal¬ 
ings of quantities integrated within a given overdensity (for 


which 

we 

drop the explicit dependence on A): 


Mgas 

oc 

Ma, 

( 3 ) 

kT 

oc 

[E(z)Ma]"/® , 


L 

oc 

^(^)2+2a/3^^+2a/3^ 



Here and throughout this work, we define L to be the in¬ 
trinsic rest-frame soft-band X-ray luminosity of the ICM 
(specifically, energies of 0.1-2.4 keV). For the kT > 5keV 
clusters in our data set, the soft-band luminosity within a 
fixed radius has a relatively mild dependence on tempera¬ 
ture, L oc {kT)°‘ with a « —0.13. This is in contrast to 
the bolometric luminosity (asoi ~ 1/2), for which the corre¬ 
sponding self-similar scaling is Lboi oc E{zy^^M‘^^} 

Where we employ log-normal probability distributions 
in this work, these are defined as Gaussian in the natural log 
of the argument, such that their standard deviations can be 
approximately thought of as a fractional scatter. 


2 DATA SET AND SPECTRAL ANALYSIS 

This work employs the same sample of 40 massive, dynami¬ 
cally relaxed clusters used to constrain cosmological models 
in Paper H. The identification of this sample is described in 
Papers I and H. Briefly, we searched the Chandra archive for 
morphologically regular clusters, as defined by a suite of im¬ 
age statistics designed specifically for the task. These mea¬ 
surements probe the sharpness of the cluster surface bright¬ 
ness peak, the alignment of a series of standard isophotes 
with one another, and the symmetry of those isophotes 
about a globally defined cluster center. In addition, a mini¬ 
mum temperature of 5 keV (excluding the core) is required.^ 
The cluster sample is summarized in Table 1. 

Our procedure for reducing and cleaning of the Chandra 
data is described in detail in Paper H. The only difference 
is that in this work we employ a more recent version of the 
Chandra analysis software and calibration files (specifically 
CIAO® 4.6.1 and CALDb"^ 4.6.2). This update, from CALDB 
4.4.10 to 4.6.2, spans a relatively important change to the 
model of the contaminant affecting the Chandra-ACIS detec¬ 
tors that, at some level, affects the analysis of observations 
as far back as 2004.® The update is, however, minor during 
the period when most of our data were observed. Comparing 
results from both versions of the analysis, we find a small 
(5 ± 5 per cent) reduction in temperatures on average, and 
negligible change in density measurements.® 

Paper II also describes in detail our procedure for 

^ Note that the 2/3 terms in the exponents of the L-Ma relation 
are identified with the slope of the kT-M/^ relation, making these 
two scaling relations less than independent where departures from 
self-similarity are concerned (see discussion by Maughan 2014). 

^ This temperature threshold is motivated by the expectation 
from simulations that the bias and scatter introduced by the as¬ 
sumption of hydrostatic equilibrium should increase for low-mass 
clusters and groups, even after selecting the most relaxed systems 
(Nagai et al. 2007). Simulations also predict a mass dependence 
of the gas mass fraction for less massive clusters, which would 
complicate the cosmological analysis of Paper II. 

® http://cxc.harvard.edu/ciao/ 
http://cxc.harvard.edu/caldb/ 

® http://cxc.cfa.harvard.edu/caldb/downloads/Release_ 
notes/CALDB_v4.5.9.html 

® Note that the cosmological results presented in Paper II are in¬ 
sensitive to an overall bias in temperature, since in that work the 
overall normalization of cluster mass profiles is calibrated using 
weak gravitational lensing data (see also Paper IV). Subsequent 
CALDB updates as of this writing (through 4.6.10a) should have 
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Table 1. Names, redshifts and masses of the relaxed clusters in our sample (see Papers I and II). Masses are from this work, using an 
updated Chandra calibration compared with Paper II. Also shown are the clean exposure time from this re-processing, as well as the 
specific Chandra observations (OBSIDs) used. 


Cluster 2 ^2500 Af25oo OBSIDs 

(kpc) (IQI^Mq) (ks) 


Abell 2029 

0.078 

648 

± 

4 

4.17 

± 

0.08 

118.9 

Abell 478 

0.088 

620 

± 

8 

3.69 

± 

0.15 

129.4 

PKS 0745-191 

0.103 

671 

± 

7 

4.74 

± 

0.14 

148.8 

RX J1524.2-3154 

0.103 

483 

± 

9 

1.77 

± 

0.10 

40.9 

Abell 2204 

0.152 

689 

± 

13 

5.41 

± 

0.30 

90.1 

RX J0439.0-I-0520 

0.208 

496 

± 

15 

2.14 

± 

0.19 

35.5 

Zwicky 2701 

0.214 

472 

± 

7 

1.85 

± 

0.08 

111.3 

RX J1504.1-0248 

0.215 

699 

± 

15 

6.02 

± 

0.38 

39.9 

Zwicky 2089 

0.235 

452 

± 

12 

1.66 

± 

0.13 

47.0 

RX J2129.6-I-0005 

0.235 

558 

± 

13 

3.13 

± 

0.22 

36.7 

RX J1459.4-1811 

0.236 

554 

± 

17 

3.07 

± 

0.29 

39.6 

Abell 1835 

0.252 

657 

± 

7 

5.21 

± 

0.16 

183.6 

Abell 3444 

0.253 

549 

± 

11 

3.04 

± 

0.19 

35.7 

MS 2137.3-2353 

0.313 

476 

± 

10 

2.11 

± 

0.13 

50.9 

MACS J0242.5-2132 

0.314 

532 

± 

31 

2.96 

± 

0.52 

7.7 

MACS J1427.6-2521 

0.318 

449 

± 

16 

1.79 

± 

0.19 

41.3 

MACS J2229.7-2755 

0.324 

471 

± 

15 

2.08 

± 

0.20 

25.8 

MACS J0947.2-I-7623 

0.345 

608 

± 

15 

4.58 

± 

0.34 

49.0 

MACS J1931.8-2634 

0.352 

573 

± 

12 

3.86 

± 

0.24 

104.0 

MACS J1115.8-I-0129 

0.355 

561 

± 

14 

3.63 

± 

0.27 

44.3 

MACS J1532.8-I-3021 

0.363 

554 

± 

11 

3.53 

± 

0.20 

102.4 

MACS J0150.3-1005 

0.363 

432 

± 

16 

1.67 

± 

0.18 

26.1 

MACS JOOll.7-1523 

0.378 

526 

± 

16 

3.07 

± 

0.29 

49.2 

MACS J1720.2-I-3536 

0.391 

529 

± 

20 

3.18 

± 

0.35 

51.7 

MACS J0429.6-0253 

0.399 

537 

± 

36 

3.35 

± 

0.67 

19.3 

MACS J0159.8-0849 

0.404 

617 

± 

18 

5.12 

± 

0.45 

62.6 

MACS J2046.0-3430 

0.423 

427 

± 

16 

1.73 

± 

0.20 

42.8 

IRAS 09104+4109 

0.442 

518 

± 

19 

3.17 

± 

0.35 

69.0 

MACS J1359.1-1929 

0.447 

463 

± 

27 

2.27 

± 

0.40 

54.2 

RX J1347.5-1145 

0.451 

809 

± 

28 

12.19 

± 

1.28 

206.5 

3C 295 

0.460 

452 

± 

22 

2.15 

± 

0.31 

90.9 

MACS J1621.3+3810 

0.461 

501 

± 

13 

2.92 

± 

0.23 

129.3 

MACS J1427.2+4407 

0.487 

477 

± 

19 

2.60 

± 

0.32 

50.8 

MACS J1423.8+2404 

0.539 

482 

± 

11 

2.86 

± 

0.19 

122.7 

SPT-CL J2331-5051 

0.576 

425 

± 

22 

2.05 

± 

0.32 

31.8 

SPT-CL J2344-4242 

0.596 

568 

± 

29 

5.01 

± 

0.78 

10.7 

SPT-CL JOOOO-5748 

0.702 

420 

± 

34 

2.29 

± 

0.56 

28.4 

SPT-CL J2043-5035 

0.723 

385 

± 

14 

1.81 

± 

0.19 

73.3 

CL J1415+3612 

1.028 

316 

± 

11 

1.44 

± 

0.15 

349.1 

3C 186 

1.063 

330 

± 

13 

1.71 

± 

0.21 

213.8 


891,4977,6101,10434,10435,10436,10437 

1669,6102,6928,6929,7217,7218,7222,7231,7232,7233,7234,7235 

2427,6103,7694,12881 

9401 

499,6104,7940 

527,9369,9761 

3195,7706,12903 

4935,5793 

7897,10463 

552,9370 

9428 

496,6880,6881,7370 

9400 

928,5250 

3266 

3279.9373 

3286.9374 
2202,7902 
3282,9382 

3275.9375 
1649,1665,14009 

11711 

3261,6105 

3280,6107,7718 

3271 

3265,6106,9376 

5816.9377 
10445 

5811.9378 

506,507,3592,13516,13999,14407 

578,2254 

3254,6109,6172,7720,9379,10785 

6112,9380,9808,11694 

1657,4195 

9333,11738 

13401 

9335 

13478 

4163,12255,12256,13118,13119 

3098,9407,9408,9774,9775 


analysing X-ray spectra for the clusters in our sample. We 
perform deprojections of the intracluster medium under two 
sets of assumptions, respectively using the nfwmass and 
PROJCT models in xspec.^ In the hrst analysis, we assume 
that the cluster mass (dark plus baryonic) is described by 
a Navarro et al. (1997, hereafter NFW) profile. The ICM 
is modeled as a series of concentric, isothermal shells, and 
is assumed to be in hydrostatic equilibrium. This analysis 
provides simultaneous constraints on the mass, temperature 
and gas density prohles of a given cluster, but has the disad¬ 
vantage of over-constraining the temperature prohle. That 
is, the shape and especially the precision of the resulting 


no effect on our analysis, given that the data employed here were 
all observed prior to 2013 (see Paper I). 

^ http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/ 


temperature profile is largely driven by the shape of the 
surface brightness profile, which has much higher signal-to- 
noise than the temperature information available intrinsi¬ 
cally in the spectra. Consequently, we use these results only 
for determining cluster mass prohles, as well as the covari¬ 
ance between mass and temperature measurements (see be¬ 
low). 

For the thermodynamic measurements which are the 
focus of this work, we instead use the “non-parametric” 
analysis described in Paper II. This analysis also assumes 
spherical symmetry, but does not require hydrostatic equi¬ 
librium. The results therefore contain no information about 
the cluster mass, but by the same token they are indepen¬ 
dent of assumptions regarding the mass prohles, and hence 
do not over-constrain the temperature prohles. Note that 
there is negligible difference between the gas density pro- 
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files recovered from the two fitting methods compared to 
the uncertainties. 

The annuli in which we extract spectra to analyse are 
chosen to provide good signal-to-noise for measurements 
of the gas density. To constrain temperatures, a more de¬ 
manding task, it is necessary to assume isothermality across 
groups of a few adjacent spherical shells. The effective reso¬ 
lution of our temperature profiles is thus lower than that of 
our density profiles from this analysis.® Figure Al in the ap¬ 
pendix shows the temperature and density profiles measured 
for each of our clusters, from both of the analyses described 
above. As already mentioned, the temperature prohles can 
be more finely binned when hydrostatic equilibrium and a 
form for the mass profile are assumed. Note that the good 
agreement between temperature profiles from the two anal¬ 
yses verifies that the NFW profile-hydrostatic equilibrium 
model is a good description of the data, within measure¬ 
ment uncertainties. 

Our primary results are based on the profiles shown 
in Figure Al. However, we also consider integrated (i.e. 
not radially resolved) thermodynamic measurements in Sec¬ 
tion 5.2. For these results, we use the same apparatus as 
the non-parametric model described above, but with some 
specializations. Specifically, this model has only two free 
temperatures, corresponding to the cluster volumes at radii 
< 0.15 rsoo and > 0.15 rsoo, with rsoo estimated from the 
NFW ht. The temperature in the outer radial bin is essen¬ 
tially identical to what we would obtain from a typical anal¬ 
ysis of the spectrum in projection, but this approach allows 
the measurement covariance between this “center-excised” 
temperature, the projected luminosity, and the spherically 
integrated gas mass to be fully captured. 

In Paper II, we identihed for each cluster a central re¬ 
gion that was excluded from the NFW model fits. By default, 
this was a circle of radius 50 kpc, although in some cases the 
excluded region is larger in order to encompass visible struc¬ 
ture in the ICM.® The motivation for excluding these data 
from the fit is to avoid biases due to local violations of hydro¬ 
static equilibrium, gas inhomogeneities, asphericity, and any 
offset between the densest gas and the center used for depro¬ 
jection (which was chosen based on the large-scale emission). 
Equilibrium is not a concern for our non-parametric fits, but 
the other issues above can potentially still bias the results 
of a spherical deprojection at small radii. In particular, we 
generically expect non-zero ellipticity or an offset of the de¬ 
projection center from the brightest pixel to result in an un¬ 
derestimate of the gas densities and/or overestimate of the 
temperatures at the smallest radii. In the following sections, 
we will present results from non-parametric fits which ex¬ 
tend all the way into the cluster centers alongside (more ro¬ 
bust) results where the central regions are discarded. Fits to 
the thermodynamic prohles will always use center-excluded 
data. 

In Paper II, we calibrated our X-ray mass measurements 
to a standard provided by the Weighing the Giants weak 
lensing analysis (von der Linden et al. 2014; Kelly et al. 


° We also leave metallicities free, although these must be binned 
even more coarsely than temperatures. 

® Point sources are always masked, regardless of their position 
within the cluster. 


2014; Applegate et al. 2014). As detailed further in Paper IV, 
this correction factor is consistent with unity (0.96 ± 0.08), 
and is consistent with being constant as a function of mass. 
For simplicity, we therefore do not apply any lensing-based 
correction to the X-ray NFW masses in this work. 


3 ADDITIONAL METHODOLOGY 

At various points in the following sections, we ht model pro¬ 
hles and trends with mass and/or redshift. This section sum¬ 
marizes the approaches we employ in these analyses. 

In Section 5.2, we constrain the multivariate scaling re¬ 
lations of X-ray luminosity, temperature and gas mass as 
a function of total mass and redshift for our sample. For 
this analysis, we approximate the measurement uncertain¬ 
ties on all quantities as log-normal, and the intrinsic scatter 
as a multi-dimensional log-normal distribution, so that the 
power law model for the scaling relations becomes a linear 
model in the transformed (logarithmic) variables. Covari¬ 
ances among the measured quantities are accounted for (see 
below), as well as covariances in the intrinsic scatter. Our 
htting method is inspired by the work of Kelly (2007), who 
proposed an efficient Gibbs sampling algorithm for htting 
a similar model with only one response variable. An asso¬ 
ciated paper (Mantz 2016) describes our implementation of 
a multivariate version of the Kelly (2007) algorithm, which 
is publically available as a package for called lrgs.^^ 
Note that generic Markov Chain Monte Carlo algorithms 
have signihcant difficulty navigating the complex parameter 
space of the elements of a covariance matrix, at least in our 
experience, making the Gibbs sampling approach of lrgs 
potentially useful for multivariate regression in general. 

The thermodynamic results that we present in this 
paper stem from the non-parametric spectral analysis de¬ 
scribed in the previous section. The uncertainties from that 
analysis directly rehect the statistical power of the data. 
However, in order to constrain scaling relations with mass 
for the entire sample of clusters, most of which do not have 
weak lensing data, we must use mass measurements from 
the hydrostatic NFW analysis. This introduces a dilemma 
without a perfect solution: we would like to use masses 
from the NFW analysis but, to avoid spuriously tight con¬ 
straints, temperatures from the non-parametric analysis. At 
the same time, it is crucial to account for the covariance 
between mass and temperature introduced by the hydro¬ 
static assumption, i.e. the fact that statistical fluctuations 
in temperatures also affect the measured masses. We deal 
with this issue for each cluster as follows. First, we iden¬ 
tify posterior samples of total mass from the NFW analy¬ 
sis, sorted into increasing order, with samples of the other 
X-ray observables from the non-parametric analysis sorted 
into increasing-temperature order. At this point, the mass- 
temperature correlation coefficient from the list of samples is 
unity. We then randomly permute elements of the mass list 
until the mass-temperature covariance in this hybrid list of 
samples matches the mass-temperature covariance internal 


http://www.r-project.org 
https://github.com/abmantz/lrgs 
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to the NFW fit (as measured directly from the NFW poste¬ 
rior samples). The result of this procedure is a set of sam¬ 
ples with marginal temperature variance matching that of 
the non-parametric fit, marginal mass variance matching the 
NFW fit, and with an orientation of the mass-temperature 
error ellipse that also matches the NFW fit. The internal 
covariance of gas mass, temperature and luminosity is unaf¬ 
fected by this process. Note that we implicitly assume neg¬ 
ligible mass-Mgas and mass-luminosity correlations at fixed 
radius; certainly we do not expect strong correlations to be 
introduced by the hydrostatic assumption apart from that 
between mass and temperature (given the use of a soft-band 
luminosity; see Section 1). 

In Section 5.1, we simultaneously fit a model for the 
mean scaled cluster density and temperature profiles, the 
intrinsic scatter between these profiles and internal to each 
as a function of radius, and the overall scaling of each profile 
with mass and redshift. The model is similar to the “Gaus¬ 
sian process” pressure profile model used by Sayers et al. 
(2013), consisting of a value for the mean profile at a set 
of scaled radii, the intrinsic covariance among profiles at 
those radii, and overall scalings with mass/redshift.^^ We 
expand on the Sayers et al. (2013) approach by marginaliz¬ 
ing over uncertainty in the mass profile of each cluster. We 
use the same approach as above to account for the mass- 
temperature covariance, but in this case the mass affects not 
just the scaling of the density and temperature, but also the 
scaling of the ordinate that the mean profile is a function of, 
i.e. A or r/r 25 oo- As above, the intrinsic and measurement 
covariances are modelled as log-normal (now a function of 
radius). In practice, we interleave Metropolis steps for the 
global scaling parameters (the exponents in the scaling re¬ 
lations) with LEGS Gibbs samples of the mean profiles and 
intrinsic covariance. 


4 THERMODYNAMIC PROFILES 

To summarize the results of our non-parametric analysis for 
the entire cluster sample, Figure 1 shows “ensemble” profiles 
of electron density, temperature, pressure (nkT) and entropy 
after applying scalings from the Kaiser (1986) 
model. Dark and light shaded regions show, at each radius, 
the 68.3 and 95.4 per cent limits of the distribution of results 
for all clusters that provide data at that radius, reflecting 
both measurement uncertainty and intrinsic scatter. These 
figures thus represent the results in a minimally processed 
form, whereas estimates of the intrinsic, cluster-to-cluster 
scatter of the profiles require us to adopt and fit a model 
(see Section 5.1). The shaded part of each profile shows the 
results when the central region of each cluster is excluded, 
while unshaded lines show the continuation of the profiles 
into the cluster centers (see Section 2). Note that the entire 
cluster sample contributes to the ensemble at overdensities 
A > 10^; however, the number of clusters where we can 
confidently measure temperatures drops quickly at larger 

Pedantically speaking, neither work employs a Gaussian pro¬ 
cess, since the models are parametrized at a finite number of 
points, between which we interpolate. Nevertheless, we perpetu¬ 
ate this incorrect terminology in the hope that future work will 
move beyond these relatively limited, discrete models. 


radii, and results for A <, 500 should be particularly treated 
with caution. 

In Figure 1, the thermodynamic profiles are plotted 
against two ordinates: the overdensity (Equation 1), and the 
physical radius in units of r 25 oo, both determined from the 
best fitting NFW profile for each cluster. The former is a 
natural choice for examining the scaling relations given in 
Equation 2, which in the self-similar model would hold at 
all A, and their intrinsic scatter. 

The following subsections describe some of the qualita¬ 
tive features of these thermodynamic profiles, and compare 
them to results in the literature. 

4.1 Density and Surface Brightness 

The self-similar assumption of a constant gas-mass fraction 
with radius leads to a simple, mass-independent scaling for 
gas density at a given overdensity, n oc E{zfiA. Figure 1 
shows that this relation approximately holds over a wide 
range in A; the dashed line in the top-left panel corresponds 
to a power law nE~^ oc A° ®®, fit to the median center- 
excised profile. Qualitatively, we see that the self-similar 
expectation is approximately met over most of the volume 
probed, in particular in the center-excised region. The de¬ 
parture from a precise A^ ® scaling can be seen more clearly 
in the /gas = Mg^s/M profiles presented in Paper II (Fig¬ 
ures 2 and 3), specifically as a shallow increase of /gas with 
radius. 

The intrinsic scatter in scaled density profiles has imme¬ 
diate consequences for mass proxies based on the soft X-ray 
luminosity, the most straightforward quantity to estimate 
from shallow data. To demonstrate, we extract (observer 
frame) 0.6-2.0keV band surface brightness profiles for each 
cluster. The expected self-similar scaling for surface bright¬ 
ness at fixed A is S oc kT K{z,T)E{zfi/{I zfi (e.g.. Pa¬ 
per I), where K is the redshift- and temperature-dependent 
K-correction. Adopting this scaling, we fit a mean profile and 
log-normal intrinsic scatter (including covariance terms) as 
a function of A. The results are shown in Figure 2. We find 
a log-normal intrinsic scatter of ~ 0.25-0.3 at intermedi¬ 
ate radii, 10^ >, A 10®, with increasing scatter at smaller 
radii. Despite the relatively large scatter, the simplicity of 
the mean model, which is well approximated by a power 
law, is appealing; in particular, these results imply that con¬ 
straints on the shape of a cluster’s mass profile can in princi¬ 
ple be obtained from surface brightness data only, by statis¬ 
tically associating isophotes with overdensities (as Figure 2 
shows, neglecting the temperature-dependent components of 
the scaling increases the intrinsic scatter only marginally).^® 
In practice, the relatively large scatter at fixed A makes this 
approach less effective than traditional X-ray proxies when 
one only wants to estimate a total mass, which is usually 
the case. 


Strictly speaking, this is true only of clusters that are morpho¬ 
logically equivalent to those in our sample, although we expect the 
surface brightness template for relaxed clusters to roughly hold for 
the general population at intermediate-to-large radii, where ear¬ 
lier studies have shown approximate self-similarity of the density 
profiles (e.g. Vikhlinin et al. 2006; Croston et al. 2008; Paper I). 
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Figure 1. Profiles of electron number density, temperature, pressure and entropy, scaled according to the Kaiser (1986) model, as a 
function of overdensity and scaled radius. The 68.3 and 95.4 per cent confidence regions at each radius are shown, where these probabilities 
encompass both measurement uncertainties and intrinsic scatter among clusters in our sample. Blue lines and shading show the results 
when excluding the central region of each cluster (as defined in Section 2), while the results shown with black lines include the centers 
(where we expect larger systematic uncertainties due asphericity). The two sets of results do not match precisely at all radii where they 
overlap because the central exclusion radius corresponds to different values of A or r/r2500 from cluster to cluster; hence, the two profiles 
differ at radii where only a subset of clusters contributes to the center-excluded results. In the upper-left and lower-right panels, the 
dashed lines are power-laws fit to the center-excised data, n oc and K oc {r/r 2 BOo)^'^^. 
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1 ^ 




Figure 2. Left: surface brightness in the 0.6-2.0keV band for each of our clusters, scaled with redshift as indicated. For clarity, 
measurement uncertainties are not shown. Right: intrinsic scatter in surface brightness at fixed overdensity. Open circles show the scatter 
when the profiles are scaled only by the redshift-dependent factor E{z)^/{1 + z)^, as in the left panel, while filled circles show the slightly 
smaller scatter achievable by including the temperature-dependent scaling terms. 


It has long been known that luminosities measured ex¬ 
cluding cluster centers (by some reasonable definition) have 
a smaller intrinsic scatter with mass than the total lumi¬ 
nosity (e.g., Fabian et al. 1994; Markevitch 1998; Maughan 
2007; Mantz et al. 2010b). For the cluster population in gen¬ 
eral this is intuitive, since the presence or absence of a cool 
cores can mean as much as factor of two difference in to¬ 
tal luminosity. The increase in surface brightness scatter at 
small radii, in part reflecting the evolution in the central 
brightness of the most relaxed clusters noted in Paper I (see 
also Santos et al. 2010; McDonald et al. 2013), shows that 
a center-excised luminosity should display reduced scatter 
with mass even within our morphologically similar cluster 
sample. As it happens, for realistic NFW concentration pa¬ 
rameters, the inner boundary of the lowest-scatter range 
noted above, A ~ 10"^, roughly corresponds to the inner ra¬ 
dius adopted by Maughan (2007) and Mantz et al. (2010b) 
when defining center-excised luminosity, 0.15 rsoo- We return 
to the subject of scaling relations using integrated X-ray ob¬ 
servables in Section 5.2. 


4.2 Pressure 



Figure 3. Scaled ensemble pressure profile as a function of 
r/rsoQ. The thick, gray line shows a generalized NFW fit to XMM 
and Planck data (Planck Collaboration 2013), while the dashed, 
red lines show the Icr intrinsic scatter about the mean profile from 
a Gaussian process fit to Bolocam data (Sayers et al. 2013). 



l/lsoo 


The shape of the ensemble pressure profile is relevant 
for reconstructing the integrated Compton-T signal from 
Sunyaev-Zel’dovich (SZ) measurements with limited reso¬ 
lution (Planck Collaboration 2014), field of view (Czakon 
et al. 2015), or sensitivity to sufficiently large scales (Mantz 
et al. 2014a). To facilitate comparison with results in the 
literature. Figure 3 shows our ensemble pressure profile as a 
function of r/rsoo. The dashed red lines in the figure show 
the la intrinsic scatter region from a Gaussian process fit 
to Bolocam data (Sayers et al. 2013; similar to the fit we 
perform in Section 5.1). The slope of the Bolocam Gaus¬ 
sian process profile is in excellent agreement with ours at 
large radii ( .> 0.5rsoo), but appears to flatten earlier at 
small radii; this difference can plausibly be explained by that 
fact that the Bolocam fit was not limited to relaxed/cool- 
core clusters. The gray, solid line in the figure is the gen¬ 


eralized NFW pressure profile for cool-core clusters from 
the Planck Collaboration (2013), which is fit to a combina¬ 
tion of XMM-Newton X-ray data at small radii and Planck 
SZ data at large radii. There is an overall offset of order 
10 per cent between our median pressure profile and the 
'KMM/Planck profile, although they agree within our sta¬ 
tistical-l-systematic uncertainties. The agreement in shape 
between our profile and the XMM/Plancfc results is qualita¬ 
tively good over the entire radial range. 

For reference, we note that our median profile is well 
described by 

( \ 
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where G is a generalized NFW function (see Na- 
gai et al. 2007) with parameters (Pq, C 500 , 7 , a,/3) = 

(1,1.19,-0.01,0.51,4.37). We stress that this parametrized 
ht should be used only as a simple description of our results, 
and only over the radial range shown in the hgure. In partic¬ 
ular, the scaled profile and its intrinsic scatter from cluster 
to cluster are modeled in greater depth in Section 5.1. 


4.3 Entropy 

The ensemble entropy profile in Figure 1 displays features 
in common with numerous other studies of clusters and el¬ 
liptical galaxies, namely a power law with radius (dashed 
line in the figure), flattening to a “floor” at small radii 
(Lloyd-Davies et al. 2000; Ponman et al. 2003; Piffaretti 
et al. 2005; Donahue et al. 2006; Morandi & Ettori 2007; 
Cavagnolo et al. 2009; Pratt et al. 2010; Werner et al. 2012, 
2014; Panagoulia et al. 2014; McDonald et al. 2014). We 
note, however, that the departure from power-law behavior 
is limited to radii which are excluded from our main results 
due to the considerations discussed in Section 2, specifically 
the possibility that the choice of center, ellipticity, or pres¬ 
ence of local density fluctuations (sloshing, cavities, etc.) 
may bias the density profiles inferred from a spherical de¬ 
projection. The ensemble profiles in Figure 1 are relatively 
conservative in this regard, excluding at least the central 
50 kpc, even when no disturbances in the ICM are visible. 

To better probe the central regions of the clusters, we 
do away with this default exclusion in this section, extend¬ 
ing the entropy profiles inwards until the radii where the 
specific features mentioned above appear to become signifi¬ 
cant. The resulting individual profiles (blue lines with error 
bars) are shown as a function of metric radius in Figure 4. 
Although the “undisturbed” region of several of our clus¬ 
ter profiles reaches as low as ~ 20kpc (~0.01-0.02r 2 oo )5 
less than half the radius of the previous default exclusion, 
there are no strong indications of departures from a power 
law. Our results do not exclude the possibility of an entropy 
floor of radius ~ 30 kpc in cool-core clusters, as generally 
agreed upon in the literature. However, they do call into 
question the validity of this concordance, given that sys- 
tematics that may affect the determination of deprojected 
density are present at these radii in even the most relaxed 
and massive clusters. We note that both Panagoulia et al. 
(2014) and Werner et al. (2012, 2014) present measured en¬ 
tropy profiles that display power-law behavior down to radii 
of <, 1 kpc ( <, 0.01 r 2 oo), respectively for samples of groups 
and highly relaxed elliptical galaxies. 

The dashed, red line in Figure 4 shows the “baseline” 
entropy power law due to purely gravitational structure for¬ 
mation found from the simulations of Voit et al. (2005), 
which is in good agreement with more recent simulations 
using several different codes (e.g. Sembolini et al. 2015). 
Here we have scaled the baseline profile according to the 
median mass and redshift of our cluster sample (see Voit 
et al. 2005), and to account for the difference between the 
cluster gas-mass fraction in the simulations and that mea¬ 
sured from our data at r 25 oo. There is good agreement be¬ 
tween the theoretical curve and the data at radii .> 500 
kpc, although the power-law slope of the observed profiles is 
somewhat shallower than the model (~ 1.15 compared with 



r (kpc) 


Figure 4. Profiles of entropy as a function of metric radius in 
our cluster sample. In this figure, the profile for each cluster is 
extended to the smallest radius where the ICM is undisturbed 
and where a spherical deprojection appears valid. With this re¬ 
striction, each of the profiles shown appears consistent with a 
power law. The dashed line corresponds to the predicted behav¬ 
ior from simple hydrodynamical simulations (without cooling or 
feedback) scaled according to the median mass, redshift and gas 
mass fraction of the observed sample. 

1.21), suggesting that we may be seeing some excess entropy 
near cluster centers. 

We emphasize that a power-law entropy profile extend¬ 
ing down to relatively small radii in this sample of dynam¬ 
ically relaxed clusters is compatible with the notion that 
cooling and feedback processes significantly affect the gas at 
small radii in groups and clusters. Indeed, our results in Sec¬ 
tions 5.2 and 5.3 indicate an important role for cooling and 
feedback within ~ 0.15 rsoo for the massive clusters studied 
here. The measured entropy profiles suggest that heating 
and cooling in the central regions of these systems proceeds 
via a continuous feedback cycle that allows the central en¬ 
tropy profiles to remain stratified, preserving the similar¬ 
ity and relatively low intrinsic scatter evident in Figure 4, 
rather than violent episodic outbursts. Within the central 
20 kpc, clear departures from spherical symmetry prevent 
robust thermodynamic profiles from being extracted from 
the data. 


5 SCALING RELATIONS 

We now turn to the scaling of the measured thermodynamic 
quantities with mass and redshift. First, in Section 5.1, we 
consider the scaling of the thermodynamic profiles, including 
their intrinsic covariance as a function of radius. Section 5.2 
then considers the scaling relations of traditional integrated 
or global cluster measurements. 

In previous work, we have stressed that scaling relations 
of the cluster population as a whole can only be reliably 
constrained when the sample selection and the underlying 
cosmological mass function are properly incorporated into 
the fitted model (Mantz et al. 2010a; Allen et al. 2011; see 
also Evrard et al. 2014; Czakon et al. 2015). This is impossi¬ 
ble in the present case, since the sample selection is complex 
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(Paper I) and depends on features which are not yet accu¬ 
rately reproduced in simulations (e.g., cool cores). Among 
the selection criteria, the kT > 5 keV cut is clearly impor¬ 
tant in the context of scaling relations. Most of our clusters 
are also initially selected from X-ray flux limited surveys, 
and must satisfy a redshift-dependent cut on central surface 
brightness, so we might expect selection effects to influence 
the density or luminosity scaling relations at some level. It 
is less clear that the selection should affect gas masses or 
center-excised luminosity measurements, which are domi¬ 
nated by the fainter parts of a cluster. While it has been 
argued that marginalizing over a flexible prior on the un¬ 
derlying distribution of masses and redshifts, as we do, may 
partially mitigate these issues (Sereno & Ettori 2015), our 
results in this section should be interpreted as strictly em¬ 
pirical, and certainly only apply to relaxed clusters. 


5.1 Scaling of Profiles 

Using the method described in Section 3, we fit a joint model 
for the mean density and temperature profiles, their in¬ 
trinsic covariance, and their scaling with mass and E{z). 
In detail, this means that at a hxed set of Na ordi¬ 
nates, {Ai}, we constrain the scaled density and tem¬ 
perature, ni[E(2)/£i(0.35)]“'^""’(M25oo/3 x 
and kTi[E{z) /£;(0.35)](M 2500 /3 x lOi"* Mq The 
slopes, P, are also free parameters of the fit. The model also 
includes a 2Na x 2Na log-normal intrinsic covariance ma¬ 
trix, encoding the marginal cluster-to-cluster scatter in den¬ 
sity and temperature at each A, the covariance of each of 
these quantities at each pair of overdensities (e.g., the covari¬ 
ance of departures from the mean density at small and large 
radii), and the density-temperature scatter covariance as a 
function of overdensity. Constraining this model requires the 
full set of 40 clusters to provide data at each overdensity. To 
ensure that this is the case, even though the mass prohle of 
each cluster (effectively, A as a function of angular radius) 
is being marginalized over, we model only the overdensity 
range 10"^'® ^ A ^ 10®'^. Table 2 summarizes the results. 
Note that it is straightforward to convert this model for the 
joint density-temperature scaling and scatter to a model for 
pressure, entropy, or some other product of powers of n and 
kT. 

Figure 5 shows the 68.3 and 95.4 per cent conhdence 
constraints on the power-law slopes of density and temper¬ 
ature with mass and E{z). Black circles show the predic¬ 
tions of the self-similar model, and dashed lines correspond 
to varying a single exponent at a time in Equation 2 (in¬ 
cluding the implicit in the expression for n). All of the 
measured slopes are consistent with the self-similar model. 

The right panel of Figure 5 shows the marginal intrin¬ 
sic scatter in density and temperature at each overdensity. 
Consistent with our analysis of the surface brightness in Sec¬ 
tion 4.1, the scatter in density is largest at small radii, and 


Note that there is no particular reason to expect departures 
from self-similar evolution to manifest themselves as a power-law 
in E{z), as opposed to some other function of redshift. However, 
this parametrization does provide a simple and convenient null 
test of whether such departures exist. 


levels off at A <, 10^. The temperature scatter shows a simi¬ 
lar trend, and is smaller than both the scatter in the density 
profile and the scatter in average temperature determined in 
Section 5.2. The latter might be expected, given that in this 
section we are taking advantage of the radially resolved, 3- 
dimensional temperature profiles. However, we caution that 
the individual temperature measurements are still effectively 
averaged over a relatively large volume (i.e. the resolution 
in radius is poor) compared to the density; this might result 
in intrinsic scatter on smaller scales, which can be probed 
in density, being suppressed in the temperature prohles. A 
comparison to realistic hydrodynamical simulations and/or 
very deep, high-resolution observations could shed light on 
this question. 

Constraints on the off-diagonal covariances of the 
model, in the form of correlation coefficients, are shown in 
Figure 6. The correlation of density at different radii appears 
to drop monotonically with their separation, approximately 
independent of their absolute position within the cluster (at 
least, within the modeled radial range). The same is true of 
temperature, albeit with larger uncertainties. This indicates 
that the overall scalings with mass and E{z) are broadly 
successful in standardizing the density and temperature pro¬ 
hles, whereas strong correlations or anticorrelations at large 
separations would suggest that perhaps the scalings should 
be radius-dependent. The density-temperature correlations 
are consistent with zero with the exception of anticorrela¬ 
tions at small separation and small radius; this presumably 
rehects local pressure equilibrium even in the presence of 
variations in the spatial extent and magnitude of cool cores. 

Overall, these results are consistent with what we would 
expect for the most relaxed clusters in the Universe, i.e. sim¬ 
ple mean scalings with mass and redshift and similar prohle 
shapes (outside the core). Expanding this kind of analysis 
to a more dynamically diverse cluster selection, most likely 
in simulations given the difficulty in measuring precise mass 
prohles for unrelaxed clusters, could prove interesting in the 
future. 

5.2 Scaling of Integrated Measnrements 

We next turn to the more traditional scaling relations of 
global or integrated quantities. For this analysis, we adopt 
rsoo, the most common choice in the literature, as the char¬ 
acteristic radius within which to make these measurements. 
Here the values of rsoo are determined by our NFW anal¬ 
ysis of each cluster, and we account for the correlation of 
mass and temperature at hxed radius imparted by the hy¬ 
drostatic assumption as described in Section 2, as well as 
the straightforward correlation in measurement uncertainty 
between mass and other quantities measured within a mass- 
dependent aperture. The thermodynamic quantities we con¬ 
sider are the spherically integrated gas mass within rsoo, 
the average temperature measured in an annulus spanning 
0.15 < r/rsoo < 1 (as determined from a 2-temperature-bin 
non-parametric deprojection; see Section 2), the soft-band 
(0.1-2.4 keV) intrinsic luminosity projected within rsoo, and 
the “center-excised” luminosity projected onto an annulus 
spanning 0.15-1 rsoo, denoted Lee. These measurements are 
plotted in Figure 7, along with a representation of their 
covariance with mass. Using the regression technique de¬ 
scribed in Section 3 and by Mantz (2016), we simultane- 
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Table 2. Constraints on the mean scaled electron density and temperature profiles of our sample, the power-law scaling indices of 
each profile with M 2500 ^-nd E{z), and the intrinsic covariance about the mean profiles. The model is given in terms of a discrete set of 
overdensities (as a proxy for radius), and the scaling is relative to 2 : = 0.35 and M 2500 = 3 x 10^^ The covariance is broken down into 
the density—density, temperature—temperature and density-temperature blocks, and expressed as marginal standard deviations along the 
diagonal and correlation coefficients off of the diagonal. 


Mean Slopes 


A 

104.6 

104.3 

104.0 

103.7 

103.4 


E{z) 

M 2500 

ln(n/cm“^) — 

4.51 ± 0.05 

-5.18 ±0.04 

-5.83 ±0.03 

-6.51 ± 0.03 

-7.18 ±0.03 n 

2.0 ±0.2 

0.03 ±0.06 

ln(fcr/keV) 2.007 ± 0.030 

2.077 ± 0.023 

2.107 ±0.014 

2.097 ± 0.014 

2.073 ±0.017 kT 

0.64 ±0.08 0.67 ±0.02 

Covariance 
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Figure 5. Left and center: constraints (68.3 and 95.4 per cent confidence) on power-law slopes of the overall scaling of gas density and 
temperature profiles with mass and E{z). The self-similar expectation is indicated by filled black circles, and dashed lines show the effect 
of changing a single exponent at a time in Equation 2. Right: marginal intrinsic scatters of density and temperature at the modelled set 
of overdensities (see Table 2). 


ously fit for the scaling of these 4 variables as a function 
of Msoo and E{z), including their 4x4 intrinsic covariance 
matrix. For completeness, we also include below some re¬ 
sults for the luminosity projected within a radius of 0.15 rsoo 
(Lcore = L — Lee), although the computationally singular 
correlations among the three luminosity measurements pre¬ 
vent us from fitting all 5 scaling relations simultaneously. 

Figure 8 and Table 3 show constraints on each of the 
scaling relation slopes. The gas mass and temperature slopes 
are all consistent with self-similarity, as expected from our 
analysis of the density and temperature profile scalings in 
Section 5.1.^® However, the scaling of total luminosity is not 

We note in particular that the steep Mgas—M slopes inferred 
from some analyses of dynamically heterogeneous systems, ex- 


consistent with self-similarity, preferring either a weaker de¬ 
pendence on redshift, a stronger dependence on mass, or 
both, consistent with a number of previous results (e.g., 
Reiprich & Bohringer 2002; Zhang et al. 2007, 2008; Mantz 
et al. 2008, 2010b; Rykoff et al. 2008; Pratt et al. 2009; 
Vikhlinin et al. 2009; Leauthaud et al. 2010; Reichert et al. 
2011; Sereno & Ettori 2015; see also the review of Giodini 
et al. 2013). This is a clearly radius-dependent phenomenon, 
with Lee being perfectly consistent with self-similarity, Lcore 
preferring stronger departures from self-similarity, and the 
total luminosity occupying a middle ground. While the un¬ 


tending down to the group scale, e.g., Pratt et al. 2009; Sun et al. 
2009) are not consistent with the results reported here for the 
most massive, dynamically relaxed systems. 
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Figure 6. Correlation coefficients of the joint density-temperature intrinsic covariance, plotted as a function of separation in overdensity 
units, n—n and kT-kT correlations appear to drop with separation in a roughly translation-invariant way. The kT-n correlations are 
all approximately consistent with zero except at small radii (large Ai A 2 ). Small horizontal offsets are introduced to make the plots 
clearer. 






Figure 7. Scatter plots summarizing the integrated thermodynamic quantities for which we fit scaling relations with M 500 and E{z). 
In each panel, the measurement covariance ellipse is shown for the most massive cluster in the sample. Shaded regions show the Icj 
predictions for a subset of the model space we explore, specifically with the power of E{z) fixed to 0.0 (for Mgas) or 2.0 (for L and Lee), 
or required to be equal to the power of M 500 (for kT). 


certainty on the evolution of I/core is substantial, milder than 
self-similar evolution is qualitatively consistent with the de¬ 
creasing central surface brightness of cool-core clusters with 
redshift seen by Santos et al. (2008, 2010) and McDonald 
et al. (2013), and in Paper 1. 

Constraints on the intrinsic covariances of the model 
appear in Figure 9 and Table 4. Our constraints on the 
marginal scatter of Mgas and kT are largely consistent with 
earlier work (see Allen et al. 2011; Giodini et al. 2013, and 


references therein); in particular the small scatter in Mgas, 
0.086 ± 0.023 is similar to the scatter found for /gas in a 
shell spanning radii 0.8-1.2r25oo in Paper II, 0.074 ± 0.023 
(see also Section 5.3; Allen et al. 2008). The kT scatter, 
0.134 ± 0.019, agrees well with results for the cluster pop¬ 
ulation at large (e.g., Arnaud et al. 2007; Vikhlinin et al. 
2009; Mantz et al. 2010b). As one would expect, given the 
intentional morphological similarity of the clusters in our 
sample, we measure a scatter in total luminosity at fixed 
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Figure 8. Constraints on the power-law slopes of the scaling relations of integrated gas mass, temperature and luminosity as a function 
of mass and E{z). The self-similar expectation is indicated by filled black circles, and dashed lines show the effect of changing a single 
exponent at a time in Equation 3. In the right panel, blue shading corresponds to the slopes for total luminosity (L), green shading 
corresponds to those for Lee, and orange shading corresponds to Leave- 


Table 3. Constraints on the power-law slopes of the scaling rela¬ 
tions of gas mass, temperature and luminosity (integrated within 
?’50o) as a function of mass and E{z). Normalizations are given 
as the natural logarithm of the listed quantities at z = 0.35 and 
A^SOO = 6 X 10^^ Mq. 



normalization 

E{z) slope 

Msoo slope 

Mgas /Mq 

31.98 ±0.02 

-0.11 ± 0.18 

1.04 ± 0.05 

fcT/keV 

2.18 ± 0.02 

0.61 ± 0.20 

0.66 ± 0.05 

L/erg 

103.70 ± 0.05 

1.20 ±0.43 

1.26 ±0.11 

Tce/ergs“i 

102.66 ±0.04 

1.82 ±0.35 

1.02 ±0.09 


mass, 0.24 ± 0.05, that is smaller than that of the popu¬ 
lation at large (~ 0.4; e.g., Vikhlinin et al. 2009; Mantz 
et al. 2010b). However, even with this selection, the scatter 
in Leave IS Still significant, 0.34lQ'Qg. The scatter in Lee is 
smaller than that in the total luminosity, although at low 
confidence, 0.17 ± 0.05. 

Although the sample employed here is relatively small, 
we are nevertheless able to place constraints on the off- 
diagonal terms relating Mgaa, kT, L and Lee- Our con¬ 
straint on the kT-L correlation coefficient is —0.06 ± 0.24. 
The M^es-kT correlation is constrained at a similar level, 
—0.18 ± 0.28, consistent with zero. In contrast, the Mgaa-T 
correlation prefers positive values, 0.431 q'3q, reflecting the 
fact that these quantities are both essentially integrals of the 
same gas density profile, but with different weighting. The 
correlation between Mgas and Lee is very strong, 0.88lo "lei 
correspondingly, the correlations of kT and L with Lee are 
similar to the Mg^s-kT and Mge.s-L correlations.^® There 
have been few previous measurements of the off-diagonal 
terms of the cluster scaling relation intrinsic scatter. Mantz 
et al. (2010b) measured the correlation of L and kT at fixed 
mass to be 0.09 ± 0.19 (not restricted to relaxed clusters), 
consistent with our findings above. Rozo et al. (2010) placed 
a lower limit on the correlation of L and M at fixed opti¬ 
cal richness (as defined for the MaxBCG catalog) of > 0.85. 

^® The strength of the A7gas—Tee correlation explains why Mantz 
et al. (2010b), using Mgas as a proxy for total mass, found a 
smaller intrinsic scatter in Lee than we do ( < 10 per cent). 


More recently, Maughan (2014) fit Lee, kT and Mgas scal¬ 
ing relations, including intrinsic scatter covariances although 
without modeling measurement covariances, to published 
data for REXCESS clusters. Those results on the intrinsic 
correlation coefficients are compatible with, and have similar 
precision to, our constraints from relaxed clusters. 

Overall, the picture that emerges is one in which the 
gas temperature and density outside of the core, roughly at 
radii 0.15-1 rsoo, behaves simply and in accordance with the 
self-similar model for these massive, relaxed clusters. In par¬ 
ticular, the strong correlation in scatter between Mgas and 
Lee requires the gas density profiles to be smooth at these 
radii (Figure 1; recall that the absence of substructure on 
these scales is a requirement of the sample selection in Pa¬ 
per I), since scatter in density is reflected disportionately in 
the luminosity compared to the gas mass. In the centers of 
relaxed clusters, we see breaking of the self-similar model, 
specifically in the form of enhanced intrinsic scatter in lu¬ 
minosity, and a preference for scaling relation slopes of core 
luminosity with mass (redshift) that are steeper (shallower) 
than the self-similar prediction. These results require the ac¬ 
tion of a non-self-similar astrophysical process such as AGN 
feedback (Fabian 2012; McNamara & Nulsen 2012, and refer¬ 
ences therein), although in general both heating and cooling 
processes may contribute to the observed trends. 

The product of gas mass and temperature, lx = 
MgevakT, has been used extensively as a mass proxy in recent 
years (e.g., Maughan 2007; Vikhlinin et al. 2009; Andersson 
et al. 2011; Menanteau et al. 2012; Benson et al. 2013; Brod- 
win et al. 2016). Given a power-law plus log-normal scatter 
model, the lx scaling relation and its scatter can be de¬ 
rived directly from our results for the Mgas and kT scaling 
relations, and their intrinsic covariance. The constraints on 
the lx slopes appear in the left panel of Figure 10, and 
are consistent with self-similarity (since Mgas and kT are, 
individually). In principle, masses estimated from such a 
combination can be no more precise than estimates from 
the individual measurements, accounting for the full intrin¬ 
sic covariance matrix.Nevertheless, we show in the right 

Note, however, that part of the motivation for using lx is 
due to the fact that, in simulations, it appears to be less sensitive 
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Correlation coefficient 


Figure 9. Marginal scatters (left) and correlation coefficients (right), respectively representing the diagonal and off-diagonal elements 
of the intrinsic covariance of integrated Mgas, kT, L and Lee at fixed M 500 and E{z). The marginal scatter of central luminosity is also 
shown in the left panel. 


Table 4. Constraints on the marginal scatters (on-diagonal entries) and correlation coefficients (off-diagonal entries), respectively 
representing the diagonal and off-diagonal elements of the intrinsic covariance of integrated Mgas, kT, L and Lee at fixed M 500 and E{z). 




kT 

L Lee 

-^gas 

0.09 ± 0.02 



kT 

-0.18 ±0.28 

0.13 ±0.02 


L 

'^•^'^-0.30 

-0.06 ±0.24 

0.24 ± 0.05 

Lee 

n 

-0.30 ±0.27 

O-SSlgag 0.17 ±0.05 


panel of Figure 10 the intrinsic scatter in mass correspond¬ 
ing to products of observables of the form where 

Y is one of Mgas, L or Lee, accounting for the intrinsic cor¬ 
relation or anticorrelation of the observables as measured 
above (e.g., Stanek et al. 2010).^® We see that choosing an 
optimal combination of observables, or simply using their 
intrinsic covariance as measured above, produces modestly 
smaller scatter compared to the individual observables in¬ 
volved. Neglecting uncertainties, the intrinsic scatter of mass 
with the combination is 0.072, compared with 

0.084 with Mgas alone or 0.21 with kT alone. Combining kT 
with Lee yields smaller gains, predictably, with the product 
^j.i.05^a95 mass scatter of 0.11, and kT^'^^ 

yielding a scatter of 0.14. While it is difficult to improve on 
the small intrinsic scatter of Mgas in the regime of our clus¬ 
ter sample, we note that the combination of temperature 
and luminosity (especially Lee) is competitive. What proxy 
or proxies are most useful in practice also depends on the 
size of the statistical and systematic uncertainties of each 


to dynamical state and closer to a power-law in mass (extend¬ 
ing to poor clusters and groups) than Mgas or kT individually 
(Kravtsov et al. 2006). Our sample, being restricted to massive, 
relaxed clusters by construction, cannot test these features. 

We have chosen to have the exponents in this combination 
sum to 2 for convenience, but note that the mass scatter is invari¬ 
ant to this choice, depending only on the relative contributions of 
kT and Y. The minimum of each curve in the figure is the scatter 
one would obtain by using the individual observables and directly 
accounting for their intrinsic covariance. 


measurement; being weighted towards smaller radii, Lee is 
generally measured with smaller uncertainty than Mgas at 
fsoo, which in turn has smaller statistical and systematic 
uncertainties than the temperature at rsoo (Figure Al). 

5.3 Scaling of /gas 

Paper II presented limited results on the scaling of /gas, 
specifically /gas in a spherical shell spanning 0.8-1.2 r 25 oo, 
and on the intrinsic scatter of /gas in a few such shells when 
assuming no scaling with mass. Here we provide a more com¬ 
prehensive treatment, differing from Paper II in the version 
of the Chandra calibration used and the consistent use of the 
same model employed in previous sections, i.e. including free 
power-law slopes with both mass and E{z) as well as intrin¬ 
sic scatter. These results provide a more complete picture 
of how the intracluster medium is redistributed, relative to 
the self-similar model, due to astrophysical processes. 

Figure 11 shows constraints on the power-law slope of 
/gas with M 2500 and the intrinsic scatter for a series of spheri¬ 
cal shells of radial width 0.4 r 25 oo- In interpreting these plots, 
it must be stressed that the results across shells are non- 
trivially correlated due to the common M 2500 values used 
in the regression, as well as the use of a parametrized mass 
profile. Nevertheless, it is clear that the gas mass fraction 
at small cluster radii increases with mass, while at larger 
radii ( .> 0.8r2soo) the slope is consistent with zero. As ob¬ 
served in Paper II, the intrinsic scatter is minimized at radii 
~ >" 2500 , clearly increasing towards smaller radii. The aver¬ 
age /gas in each shell closely traces the differential profiles 
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Figure 10. Left: Constraints on the power-law slopes of Yx = MgasfcT as a function of mass and E(z), as in Figure 8 . Right: intrinsic 
scatter in mass from products of observables of the form kT°‘Y^~°‘, where Y is one of Mgas, L or Lee. 


shown in Paper II (Figure 2), while the power-law slope 
with E{z) is consistent with zero at the ~ la level in each 
shell. Table 5 shows the constraints on the fit parameters 
for a subset of the shells, as well as for /gas in full spheres 
bounded by r 25 oo and rsoo (the latter from Section 5.2) fit 
with the same model. 

The fact that the gas mass fractions of clusters increase 
with radius and are well below the cosmic average in their 
inner regions (e.g. Allen et al. 2002, 2004, 2008; Vikhlinin 
et al. 2006; Paper II) provides some of the clearest evidence 
for feedback (even more so in galaxies; e.g. McNamara & 
Nulsen 2012). In contrast, hydrodynamic simulations with¬ 
out cooling or feedback predict a constant gas mass fraction 
as a function of mass and radius (Borgani & Kravtsov 2011 
and references therein). Beyond this simple comparison, the 
results in Figure 11 and Table 5 complement those of Sec¬ 
tion 5.2 in demonstrating that the gas distribution at small 
cluster radii has a larger scatter and greater dependence 
on mass than the gas properties at larger radii (in this case 
gas mass rather than luminosity). Interestingly, even though 
/gas in the shell radially centered at r 25 oo is constant with 
respect to M 2500 , /gas integrated in a sphere bounded by 
^2500 displays an increasing slope with mass, 0.10 ± 0.05, 
due to the influence of the central regions. When integrat¬ 
ing over the larger region bounded by rsoo, this slope is 
consistent with zero. Note that we do not expect this re¬ 
sult to apply to significantly less massive clusters or at the 
group scale, where the astrophysical processes that break 
self-similarity in the centers of our clusters may show their 
influence at larger radii (e.g. Sun et al. 2009). We also note 
that both of the spherically integrated measurements retain 
a larger intrinsic scatter, due to the inclusion of cluster cen¬ 
ters, than the asymptotic value of ~ 0.07 seen in the shells at 
r .> 0.6r25oo- As noted in Paper II, the small intrinsic scat¬ 
ter in the 0 . 8 - 1 . 2 r 25 oo shell, combined with the typically 
good signal-to-noise at these radii in Chandra data, directly 
translates into improved cosmological constraints compared 
with a larger-scatter observable. The independence of mass 
of the /gas in this shell is appealingly simple, although in 
principle a mass dependence can easily be incorporated into 
the full cosmological model (see discussion in Paper II). 


6 CONCLUSION 

We present profiles of density and temperature (and the de¬ 
rived quantities pressure and entropy) for a sample of 40 
massive, dynamically relaxed galaxy clusters, and provide 
the first constraints on the scaling with mass and redshift of 
these profiles, including their intrinsic covariance as a func¬ 
tion of radius. We also fit scaling relations of traditional, 
integrated quantities: gas mass, average temperature, total 
soft-band luminosity, and center-excluded luminosity. This 
analysis includes the full multivariate intrinsic covariance 
matrix, and our results represent the first constraints on 
some of these cross terms. 

Whether in terms of profiles or aperture-integrated 
measurements, our results support a picture in which the 
ICM in these massive, relaxed clusters follows self-similar 
scaling laws, with the exception of their innermost centers 
(radii <, 0.15 rsoo). The luminosity (i.e., gas density) in these 
central regions evolves less strongly, and has a stronger de¬ 
pendence on mass, than the self-similar prediction, consis¬ 
tent with the physical effects expected from the development 
of cool cores and heating by central AGN. Even in this mor¬ 
phologically selected sample of relaxed clusters, cluster cen¬ 
ters dominate the intrinsic scatter in X-ray luminosity at 
fixed mass and redshift. The center-excised luminosity has 
a smaller scatter, ~ 15 per cent, comparable to that of the 
average temperature. For these clusters, gas mass has the 
smallest intrinsic scatter, ~ 8 per cent, although an optimal 
combination of gas mass and temperature can also serve well 
as a total mass proxy. Consistent with the results above and 
in Paper II, we find that a shell spanning ~ 0.8-1.2r25oo is 
near optimal for cosmological studies using /gas, due to the 
small scatter of the gas mass fraction at these radii. 

Within the centers of our clusters, the measured 
trends in luminosity and the gas mass fraction indicate re¬ 
distribution of the ICM by non-gravitational processes such 
as cooling and feedback. At the same time, the measured 
entropy profiles remain regular, with small scatter, down to 
the smallest radii where deprojected measurements can be 
robustly made (20 kpc ~ 0.01 r 2 oo)- The results suggest that 
the heating and cooling processes at work in the centers of 
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Figure 11. Constraints on the power-law scaling with mass (left) and intrinsic scatter (right) of /gas measured in spherical shells. 
Horizontal bars indicate the radial extent of the regions where /gas is measured, while vertical bars show 68.3 per cent confidence 
intervals. Note that the plotted points are non-trivially correlated; nevertheless, the gross trends are clear. As a function of increasing 
radius, out to ~ 1.2r2500i the slope of /gas with mass systematically decreases, becoming consistent with zero at r > 0.8r2500- The 
intrinsic scatter decreases to a minimum at ~ ^25001 appears to flatten beyond that. 


Table 5. Constraints on the gas-mass fraction and its scaling properties in various spherical shells. The first two columns show the inner 
and outer radii of each shell in units of r2500 (however the outer radius for the last line is actually rsoo rather than a strict multiple of 
^250o)- The third column shows the model constraint on /gas in each shell at z = 0.35 and M 2500 = 3 X 10^^ Mq, and columns 4-5 show 
the constraints on power-law indices with mass and E{z). Column 6 shows the intrinsic scatter in ln(/gas) at fixed mass and redshift. 


Tin rout /gas E{z) slope M 2500 slope scatter 


0.0 

0.4 

0.082 

± 

0.002 

-0.38 

± 

0.27 

0.17 

± 

0.07 

0.18 

± 

0.02 

0.4 

0.8 

0.105 

± 

0.002 

-0.14 

± 

0.17 

0.09 

± 

0.05 

0.09 

± 

0.02 

0.8 

1.2 

0.129 

± 

0.002 

-0.15 

± 

0.17 

0.00 

± 

0.05 

0.06 

± 

0.02 

0.0 

1.0 

0.100 

± 

0.002 

-0.17 

± 

0.18 

0.10 

± 

0.05 

0.11 

± 

0.02 

0.0 

~ 2.2 

0.129 

± 

0.003 

-0.11 

± 

0.18 

0.04 

± 

0.05 

0.09 

± 

0.02 


massive clusters are tightly coupled, allowing the gas to re¬ 
main stratified on these scales. 
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APPENDIX A: INDIVIDUAL CLUSTER 
PROFILES 

Figure Al shows the deprojected temperature and elec¬ 
tron density profiles of each individual cluster in our sam¬ 
ple, as determined from the non-parametric spectral fits 
(red/orange crosses in each panel) and the fits that assume 
an NFW mass profile (blue shaded boxes and lines with er¬ 
rors). Note that the choice of a deprojection center to opti¬ 
mize symmetry of the emission on large scales may result in 
biases in the density and temperature at very small radii. For 
this reason, the cluster centers are excluded from the NFW 
fits, and there are no corresponding solutions for densities 
and temperatures in those regions (from the NFW model). 
Confidence intervals at the 68.3 and 95.4 per cent level are 
shown, with the exception of the density profiles from the 
NFW model, where only 95.4 per cent intervals are shown. 
Vertical dashed lines show the best-fitting value of r 25 oo for 
each cluster. Radii expressed in kpc and densities are de¬ 
pendent on our adopted reference cosmological model. See 
Section 2 for more details. 
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Figure Al. 
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Figure Al - continued 
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Figure A1 - continued 
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Figure Al - continued 
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